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ABSTRACT 

Electron spins in quantum dots under coherent control exhibit a number of novel feedback processes. Here, we 
present experimental and theoretical evidence of a feedback process between nuclear spins and a single electron 
spin in a single charged InAs quantum dot, controlled by the coherently modified probability of exciting a trion 
state. We present a mathematical model describing competition between optical nuclear pumping and nuclear 
spin-diffusion inside the quantum dot. The model correctly postdicts the observation of a hysteretic sawtooth 
pattern in the free-induction-decay of the single electron spin, hysteresis while scanning a narrowband laser 
through the quantum dot's optical resonance frequency, and non-sinusoidal fringes in the spin echo. Both the 
coherent electron-spin rotations, implemented with off-resonant ultrafast laser pulses, and the resonant narrow- 
band optical pumping for spin initialization interspersed between ultrafast pulses, play a role in the observed 
behavior. This effect allows dynamic tuning of the electron Larmor frequency to a value determined by the pulse 
timing, potentially allowing more complex coherent control operations. 
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1. INTRODUCTION 

The control of a single electron spin in a semiconductor quantum dot (qd) using single, ultrafast optical pulses 
has emerged as a promising path towards the development of high-speed, optically driven quantum information 
processors. One remaining obstacle for using QD spins as qubits is the need to compensate for the drift of 
random nuclear fields. This problem may be alleviated by one of several recently discovered nuclear feedback 
phenomena observed in QDs under coherent control. Such phenomena occur in electrically controlled double 
QDs, in which transition processes between electron singlet and triplet states allow the manipulation of interdot 
nuclear spin polarizations, improving coherent controlP^ In single QDs under microwave control, nuclear effects 
dynamically tune the electron spin resonance to the applied microwave frequencyPl Tuning effects are also 
observed in two-color continuous-wave (cw) laser experiments in self-assembled QDs, in which the appearance of 
coherent electronic effects such as population trapping are modified by nonlinear feedback processes with nuclear 
spinsJ3H Finally, nuclear spins have been shown to dynamically bring ensembles of inhomogeneous QDs into 
spin-resonance with a train of ultrafast pulses 00 

Recently, similar nuclear feedback phenomena have also been observed for a single self-assembled QD under 
pulsed ultrafast controlP The effect is evident when measuring the familiar "free-induction decay" (fid), equiv- 
alent in this context to a Ramsey interference experiment, of a single spin in a single QD under pulsed control. 
The Larmor frequency of the electron spin is dynamically altered by the hyperfine interaction with QD nuclei; 
the nuclear polarization is in turn altered by the measurement results of the FID experiment. The result is a 
feedback loop in which the nuclear hyperfine field stabilizes to a value determined by the timing of the pulse 
sequence. The experimental observation is treated in more detail in Ref. 9. In what follows, we elaborate on a 
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Figure 1. (a) Sample geometry; a single QD is embedded between AlAs/GaAs mirrors; a transverse magnetic field 
B is applied, (b) Level structure and optical interactions used in the experiment. Q p (green) labels the CW optical 
pump, Sleff (red) labels the pulsed coherent rotation. Wavy lines indicate spontaneous emission; the photons collected for 
measurement are indicated by the dashed box. 

mathematical model to explain the effect. We also show how this model can be used to describe related hysteretic 
optical pumping phenomena and non-sinusoidal spin-echo fringes. 

These observations and this model could have strong impact on the viability of optically controlled QDs for 
quantum information processing. Prior to the discovery of nuclear control effects, nuclear spins were generally 
regarded as a problem for quantum information processors, since they cause adverse effects such as inhomogeneous 
broadening and non-Markovian decoherence processes. While useful roles such as long-lived quantum memory 
were envisioned^ their practicality was questionable due to limited amounts of nuclear polarization and memory 
loss due to nuclear spin diffusion. The present results, along with similar workp^ show that nuclear spins may 
provide novel methods for the dynamic tuning and locking of electron spin resonances for electrons trapped in 
QDs. This may benefit the stability of future coherent devices based on QDs. 

In Sec. [2] we briefly describe the experimental circumstances and, and in Sec.[3j we show how nuclear feedback 
appeared in observations. Section [4] then describes the physical processes which we believe are most likely to 
account for the observations. With these assumed, we present in Sec. [5] a simple mathematical model that 
replicates the experimental observations. We then conclude with some prospects for how this work may impact 
quantum information processing. 

2. EXPERIMENT 

The experimental apparatus is similar to that described in Refs.[ll]and[T2j we refer the reader to these references 
for more details of the experiment, including apparatus schematics. 

2.1 Quantum Dot Sample 

The sample is grown by molecular beam epitaxy to contain roughly 2 x 10 9 cm -2 self-assembled InAs QDs, 
situated roughly 10 nm above a (5-doping layer of ~ 4 x 10 9 enr 2 Si donors. About one-third of the QDs were 
negatively charged by the 5-doping layer, and a single charged QD was spectroscopically isolated at around 
940 nm for the experiments. As shown in Fig. [lja), the QDs were embedded at the center of a planar GaAs 
microcavity (Q ~ 200) with 24 and 5 pairs of AlAs/GaAs A/4 layers in the bottom and top mirrors, respectively, 
which served to direct the QD emission towards the collection lens and decrease the optical power required for 
spin rotation. As shown in Fig. [ljb), the two spin states of the QD electron, \\) and |f), are split with a Larmor 
frequency S c = g c ^B c ^t/h. Except where noted, all data presented here uses a B cxt = 4 T magnetic field, which 
leads to S c /2ir = 25.3 GHz, The field is applied in Voigt geometry, perpendicular to the optical axis and sample 
growth direction; see Fig. [lja). In this geometry, the electron-spin ground states each couple optically to both 
heavy-hole trion states litj-D-) an d I41MT)- These trion states each contain two electrons in a spin-singlet and an 
unpaired heavy hole. 



2.2 Ultrafast Coherent Control 



At the core of our experiment is the coherent control of a single electron spin using ultrafast pulses. We coherently 
manipulate the electron spin state by applying to the QD a circularly-polarized, 4-ps optical pulse that is detuned 
by ~ 150 GHz below the exciton transitions. The pulse rotates the spin approximately around the optical axis, 
allowing superpositions of \\.) and |"f)PH2l 

There are two equivalent points of view as to how the pulse accomplishes a coherent rotation. One point 
of view, described in Refs. [T2j[l3] and elsewhere, is that the pulses work via stimulated Raman transitions. 
Physically, a train of pulses from a modelockcd laser may be considered to have a comb of frequencies; two peaks 
of such a comb can provide the two elements of a Raman transition. The train of pulses (in the time domain) 
which leads to the frequency comb (in the frequency domain) will, if their timing is tuned to the Larmor period of 
the single electron spin, always rotate the electron spin with the same phase, constructively achieving a rotation. 
If the pulse power or oscillator strength is high enough, as it is in our case, substantial rotations can occur with 
just a single pulse, and the relative timing between just a pair of pulses can achieve complete qubit control.^ 
Mathematically, the action of the pulse may be derived using adiabatic elimination of the excited states, in which 
time-derivatives of excited states in either the Schrodinger or master equation are set to zero. As a result, the 
single-pulse rotation angle is approximated as 

(i) 

where A is the detuning and £7(i) is the time-varying optical Rabi frequency, proportional to the pulse's electric 
field. 

The mathematical procedure of adiabatic elimination includes two approximations: the adiabatic approx- 
imation that, due to the large detuning, the system remains in a "slowly" changing ground state; and the 
approximation inherent in first-order perturbation theory. Only the first approximation is used in the point of 



view for ultrafast coherent control as an AC-stark shift, indicated in Ref. 14 and elsewhere. For this picture, we 
note that the pulse length is much shorter than the Larmor precession period. On this timescale, the electron 
ground state is effectively degenerate, and as a result there exists an optically dark ground state (for example, 
|t) — |4), neglecting normalization), and an orthogonal optically bright state (for example, |f) + \i))- The de- 
tails of these states will depend on the selection rules of the dot, which in turn depends on details of the dot's 
symmetry. In this point of view, the arrival of the pulse "slowly" mixes the bright state and the excited state to 
which it is optically connected. Because of the large detuning, the adiabatically changing mixed state reduces 
in energy by 



A 2 



(2) 



This is the AC-stark shift of the bright state. A spin beginning in a superposition of bright and dark states 
will therefore see its bright state component shifted by the time integral of <5w(i), thereby accomplishing a spin 
rotation. Note that the first-order approximation of / 8oj(t)dt in the limit A ^ fl(t) provides the rotation angle 
obtained from the adiabatic elimination point of view, Eq. ([I]) . Also note that the nonlinear dependence on the 
rotation angle 9 with pulse power, as reported in Ref. [l2j may be described by the time- integral of Eq. ([2]). 

2.3 Pulse Sequences 

Before being manipulated by a sequence of rotation pulses, the spin is first initialized into the ground state |f) by 
optical pumping [Fig.jljb)]. A pulse of duration T p = 26 ns from a narrow-band laser, gated with an electro-optic 
modulator (eom), drives the \l) <-> lit; ft) transition at rate O p . Spontaneous decay then shelves the spin into 
|t) within a nanosecond timescale. The optical pumping step also serves to measure the spin-state resulting 
from the previous rotation sequence: if the spin was rotated to \l) then a single photon will be emitted from the 
14-titl") ~~ ^ It) transition as the spin is re-initialized. This photon is spectrally filtered by a double-monochromator 
and detected using a single-photon counter. 

In this work we describe the results of three different pulse sequences, summarized in Fig. [2j 
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Figure 2. The three pulse sequences employed, (a) One-pulse, optical pumping sequence. After the optical pump (green 
pulse) polarizes the spin to final polarization Sf , as shown by the black line, an ultrafast it pulse (red) inverts the 
polarization, and the sequence is repeated. The polarization Sf depends on the wavelength A of the pumping laser, which 
allows an absorption spectrum, (b) Two-pulse, fid or Ramsey interference sequence. The optical pump (green pulse) 
polarizes the spin from Sf to a final polarization Sf, after which it is rotated to the Bloch-sphere equator by an ultrafast 
7r/2 pulse (red). A second n/2 pulse a time r later rotates the spin to a location on the Bloch sphere with ^-projection 
Si that depends on r, and the sequence is repeated. During optical pumping, Sf — Sf is measured, (c) Three-pulse, 
spin-echo sequence. The sequence is the same as (b), except an optical 7r-pulse is inserted approximately halfway between 
the 7r/2 pulses, refocusing nuclear Overhauser shifts in the electron spin's precession. The z-projection of the spin after 
the last pulse, Sf , depends on the offset r of the n pulse from the center of the 7r/2 pulses. 



2.3.1 One-pulse Sequence: Optical Pumping 

In a one-pulse or "optical pumping" sequence, illustrated in Fig. [2^a), a single ultrafast it pulse is applied; this 
pulse is also gated by an eom. This sequence is designed to give maximal signal during optical pumping, allowing 
measurements such as the spectral absorption while scanning the pumping laser across the resonance. 

2.3.2 Two-pulse Sequence: Ramsey interferometry/Pree Induction Decay 

In a two pulse sequence, illustrated in Fig.J^b), the spin is first initialized into After the first 7r/2 rotation, 
the electron spin freely precesses around the equator of the Bloch sphere at its Larmor precession frequency, 
including any Overhauser shifts. A second 7r/2 rotation introduces a projection onto the z-axis which varies 
sinusoidally with r; that projection is measured during the subsequent initialization step as the entire sequence 
is repeated. After averaging, the resulting trion count rate is expected to show sinusoidal fringes as a function 
of r. 

In practice, the two timed pulses are achieved by first splitting optical pulse with an unbalanced Mach-Zender 
interferometer and then gating single pulses with an EOM. A retroreflector on a computer-controlled stage allows 
the time delay r between the two rotation pulses to be scanned by up to 300 ps. 
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Figure 3. Hysteretic absorption trace observed in a single-pulse experiment, while scanning the wavelength A of the 
pumping laser down and then up. 

2.3.3 Three-pulse Sequence: Spin-Echo 

In the spin-echo experiment, illustrated in Fig. [2^c), the three separate pulses from the mode-locked laser are 
each gated by multiple eoms on two split optical paths. The two paths have different optical attenuations, 
assuring a7r/2 — it — tt/2 Hahn echo sequence. We denote the time between the first 7r/2 rotation and the ir 
rotation T, which is necessarily an integer multiple of the pulse period of the mode-locked laser, 13 ns. However, 
the optical path for the middle 7r-pulse is sent through a separate optical path with variable length tuned by the 
computer-controlled stage and retroreflector, allowing 7r/2-to-7r-pulse separation to be modulated toT-r while 
the 7r-to-7r/2-pulse separation is T + r. Fringes are therefore observed as a function of r at twice the Larmor 
frequency. The decay of these fringes as a function of T allows a measurement of T 2l described elsewhere.^ 
Here, the focus is on the behavior with respect to r. 

3. OBSERVATIONS 

Nuclear feedback effects are observed during the application of all three sequences, but they manifest in different 
ways during each sequence, as we now describe. 

3.1 Optical Pumping 

If observing the count rate as the wavelength of the optical pumping laser is scanned, one expects, in the absence 
of nuclear effects, a symmetric absorption lineshape. This is not observed, however; instead, an asymmetric 
lineshape is observed. Most notably, this lineshape shows hysteretic behavior, changing its shape depending on 
the direction of the scan and the history of the experiment. This behavior is shown in Fig. [3} Such effects were 
also observed in experiments is which the repumping was achieved with a second CW-laserpHU rather than an 
ultrafast pulse as in our case. Either way, a nuclear origin to the hysteresis is evident from the long timescale at 
which the hysteretic features survive. 

3.2 Free Induction Decay 

For the FID experiment in the absence of any electron-nuclear spin feedback mechanisms, the nuclear spins 
would be expected to fluctuate randomly on a timescale slow compared to the Larmor precession, leading to 
random Overhauser shifts of the electron's Larmor frequency due to the contact-hyperfine interaction. Repeated 
measurements of the same spin with different Larmor frequencies would lead to the sinusoidal fid fringes decaying 
with a Gaussian shape on a timescale T 2 *, known to be on the order of nanoseconds from theoretical calculations. 
However, such a Gaussian decay was not observed. Figure [4] shows the result of the FID experiment. The top 
three traces show the fringes seen as the delay r is increased, and the bottom three correspond to decreasing r. 
The oscillatory fringes, rather than decaying, evolve into a sawtooth pattern at high values of r. Similar data 
were observed elsewhere as well: see the top trace of Fig. 3b of Ref. [15] This pattern also shows hysteresis 
depending on the direction in which r is scanned. Fig. mb) illustrates the result of switching the scan direction. 




3.2.1 Spin echo 

The principal goal of the spin-echo measurement is to observe fringe visibility vanish as T is increased, allowing 
a measurement of T2. Such results are reported in Ref. [II] Not reported there are notable effects observed 
in the shape of the echo fringes as a function of r, which we qualitatively explain in the following sections. A 
qualitative model will have to suffice, since the effects themselves depend on the sample, the sample history, 
and the averaging time, and hence can only be qualitatively reproduced. Examples are shown in Fig. |5ja,b). 
The two principal effects are non-sinusoidal "spikiness" of the echo fringes as a function of r [Fig. |5ja)], and a 
decay of the fringe count rate to a constant value substantially higher than the average count-rate [Fig. [5^b)]. 
Neither effect can be explained from a simple inhomogeneous broadening model, and hence, in ensuing sections, 
we provide an explanation for these effects in terms of nuclear feedback effects. 



4. PHYSICAL PROCESSES DRIVING NUCLEAR POLARIZATION 

These data result from two competing processes, which "push" and "pull" the nuclear polarization and therefore 
the hyperfine (or Overhauser) shift in contrary directions. The "push" process is due to trion excitation or decay, 
while the "pull" process is due to nuclear spin diffusion. In what follows, we discuss the physical origins of these 
processes and introduce a simple differential equation whose stable solutions model our data well. 

Let us begin by discussing the single-electron ground state of the QD, in which the electron spin, effectively 
a dipole, sees primarily a simple Zeeman splitting 

#c-z = g c ^B S z = h6 c S z . (3) 

Note that, in our convention, the (large) magnetic field is in the z direction. The growth (and predominant 
uniaxial strain) direction is taken to be the x direction. For the electron spin, this strain plays little role, since 
a pure 5 = 1/2 dipole does not couple to the electric field gradients introduced by strain. 

The largest magnitude hyperfine interaction in a QD is the Fermi contact hyperfine interaction, 

H cM = Y,A n S-lm (4) 

n 

where S is the spin operator for the electron (S — 1/2) and I is the spin operator for the nth nucleus (/ = 3/2 
for all isotopes in GaAs, and 1 = 9/2 for the two isotopes of In.) Here the coupling constant A n varies across 
the dot and by isotope; the sum A = J2 n A n is on the order of 100 /ieV. The portion of this Hamiltonian 
proportional to S z is the dominant source of hyperfine shifts in a QD; we denote this shift for some particular 
nuclear configuration as f2 = J2 n A n In/h- 





Figure 5. Experimental results for echo fringes and theoretical model, (a) At low times r, "spiky" fringes are observed. 
These have clear deviation from sinusoidal behavior, with smoother curvature on the top of the oscillations and sharp 
cusps on the bottom. This particular data was taken at magnetic field B = 3 T with a 7r-pulse delay of T = 78 ns. (b) At 
longer times t, the fringes are often observed to decay to a final count-rate which is not the expected "middle" of the 
oscillations, but rather the maximum count-rate. This particular data was taken at magnetic field B = 4T with a 7r-pulse 
delay of T = 13 ns. (c) Mathematical model explaining both spikey data and decay to a high value. The top panel shows 
the t — > oo solution to Eq. (161, i.e. Eq. (171, in the colorscale. Superimposed in white is the average polarization (fi) for 
each r. The resulting average countrate (C3(Q)) is shown in the lower panel. The timescale and Larmor frequency have 
been chosen not to match the data, but rather to easily show spiky oscillations and a decay to a high count-rate in a single 



plot. This particular model uses D/k - 
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Since nuclear spin-flip rates during trion decay are critical for our model, we must examine not the overall 
strength of the hyperfine interaction, but rather the probability of it nipping a nuclear spin during trion decay. 
The lowest order nuclear-spin-flip rate occurs in second order perturbation theory, which yields an estimate that, 
during trion decay, the probability of a contact-hyperfine- induced nuclear spin flip at nucleus n is \A n / g[i B \ 2 , 
where g[i B is the electron Zeeman energy which must be compensated for in the trion recombination process. 
The Zeeman energy is on the order of 30 GHz x h, putting this probability at roughly 10~ 8 to 10 -9 per nucleus, 
or with a timescale of about 100 ms to 1 sec. This process may play some role during trion emission, but to 
answer whether it is the dominant transition or not, we must look at the other hyperfine interactions as well. 

Before discussing the hole, we note that the QD ground state also features the electron-nuclear dipole-dipolc 
interaction, 

H V-^o. / 3(r».S)(r».I tt )-r ? 2 t (S-I») \ 

n * n I 

where r„ = r — R„ for (assumed fixed) nth nuclear position R„ and the brackets (•) now refer to an average 
over the electron wavefunction. The integral over the QD wavefunctions make this term very small for the s- 
like symmetry of the ground state Bloch wavefunction of a conduction electron close to the T-point in GaAs. 
Certainly this interaction will be smaller than the contact term for flip-flop terms such as S~ I + , and for these 



terms dipolar processes are clearly not competitive with the Fermi contact term. This Hamiltonian does admit 
terms of the form S Z I X , which in-principle conserve the Zeeman energy. Still, however, these terms are expected 
to be many orders of magnitude smaller than the contact hyperfine term for s-shell conduction electrons. The 
S Z I X term, for example, is proportional to 

/ sin 26 cos 4> 
\ 7* 

which vanishes for a pure hydrogenic s-shell wavefunction and nearly vanishes for a GaAs conduction electron. 

The heavy hole of a trion is a different situation, however, for two reasons. The first reason is that the strain- 
split p-like configuration of the valence holes allows averages like Eq. Q to survive. The average of 1 /r 3 , where 
r is the distance from a nucleus to a hole, is a nontrivial number to calculate from first principles, especially in 
the presence of strain. For a back-of-the-envelope estimate, (1/r 3 ) should be estimated as the size of a Wannier 
function on a single atomic site for the extended hole wavefunction. More detailed calculations, such as that 
in Ref. |16[ put the strength of the heavy-hole nuclear dipole coupling in a GaAs quantum well or QD, summed 
over all nuclei, in the vicinity of 10 /zeV. The average effective field of a heavy hole at an arsenic nucleus, for 
example, is therefore on the order of this number divided by Nhr/ n , where N ~ 10 6 is the number of nuclei and 
for Arsenic 7„«5x 10 7 (T sec) -1 . Using these numbers, the local magnetic field seen by a nucleus due to the 
dipole moment of the hole is on the order of 10 G. We denote this field B^K 

The second important reason the hyperfine interactions with the hole are different from the electrons is that 
light-hole/heavy- hole mixing plays a substantial role. In particular, in addition to the Zeeman term, there is a 
strain Hamiltonian, combining to form a Hamiltonian of the form 

H h = aJ 2 x +bJ z , (7) 

for J = 3/2. The coefficient a is proportional to the amount of in-grown strain in the ^-direction and b is 
proportional to the magnetic field. In our sample, a > d. Correspondingly, what we generally refer to as a 
heavy-hole trion is, at these fields, a combination of a heavy and light hole, and its dominant spin direction in 
the ground state is at some angle between the x and z, closer to the a;-direction. Suppose this direction is labeled 
fo; then the dipolar coupling between the hole and nucleus n contains the hole-Zeeman-energy-conserving term 
(J • ro)I x , which will in general not vanish in this low symmetry situation. In fact, as we discuss below, this 
term is only a bit more than an order of magnitude weaker than the Fermi contact term. Crucially, though, 
the (J • ro)I x term is nearly energy- conserving, and therefore in second-order perturbation theory, the energy 
denominator is very small, on the order of the nuclear Zeeman energy of a few MHz instead of the electron 
Zeeman energy of a few GHz. As a result, the overall transition rate during trion emission is at least an order of 
magnitude larger than that due to the contact hyperfine interaction. Fermi's golden rule in second order allows 
an estimate of the rate T n at which a trion hole nearly polarized along the sample growth axis (orthogonal to 
the magnetic field) randomly flips the nth nuclear spin during spontaneous emission, with the nuclear Zeeman 
energy compensated by the broad width of the emitted photon (7/2-77 ~ 0.1 GHz) and the photonic density of 
states negligibly changed by the Zeeman energy of the nucleus. The result is r„ w (B^ / B ) 2 j sa 1/(10 ms). 

Although we conclude that hole-hyperfine interactions are responsible for the nuclear polarization effects in 
our experiment, this conclusion should not be construed to mean that electron hyperfine interactions play no 
role at all, or never affect nuclear spin dynamics in QDs. They do and certainly can, especially in experiments 
in the Faraday geometry where the strain and field directions are parallel, and the hole dipolar interaction plays 
a lesser role. However, we believe that the particular phenomena we report here are most likely explained by 
the hole dipolar interaction. The conclusion that hole dipolar interactions contribute predominantly to nuclear 
feedback effects in optically controlled InGaAs QDs was also reached in Ref. [5j more detailed discussion appears 
in Ref. El 

The "pull" process of our model is the presence of nuclear spin diffusion, a process known to limit nuclear 
polarization rates in QDsP^ By "nuclear spin diffusion," we are referring to the motion of nuclear spin polarization 
due to the nuclear dipole-dipole interaction. Technically speaking, diffusion in the most classical sense is not a 
complete description of the loss of nuclear polarization in a QD. The reason is that the hyperfine field of the 
electron introduces a magnetic field gradient, which makes dipole-dipole induced changes in the nuclear field 




energy-nonconserving. Polarization may still move and change, however, due to the finite dipolar width in the 
large three-dimensional bath of spiniP^ anc i due to fluctuations of the gradient.^ In general, the magnitude of 
the nuclear dipole-dipole coupling sets a spin-diffusion timescale on the order of milliseconds. However, the more 
complex processes occurring in the presence of the hyperfine field gradient allow for a hierarchy of timescales. 
The fastest, millisecond timescales would describe nuclei approaching a quasi-equilibrium dipolar configuration, 
which may depend on the history of magnetic interactions inside the dot. We denote the timescale to reach 
this quasi-equilibrium Eventually, diffusion should cause nuclear polarization to "leak out" of the QD and 
reach thermal equilibrium. This timescale may be far longer, however, on the order the nuclear T\ rate, which 
is observed to be on the order of 20 minutes in our system. 

We therefore see that our experiment features five very distinct timescales, and our model will depend 
crucially on the assumption that they remain distinct. The fastest timescale is the dynamics during picosecond 
pulses, which play no role for the nuclei; hence the electron rotations are considered instantaneous. The next 
slower timescale is that of the pulse sequence and resulting electron-spin dynamics, repeated continuously with 
a repetition period on the order of hundreds of nanoseconds in the one- and two-pulse experiments, and on 
the order of 2T in the spin-echo experiments. This is still much faster than the third timescale, that of the 
millisecond-speed nuclear dynamics. Next, the averaging timescale of the measurement is much longer, on the 
order of several seconds, allowing the nuclei substantial time to reach quasi-equilibrium. The longest timescale, 
that of nuclear T\, is far longer still, and plays no role in our model. 

5. MATHEMATICAL MODEL 

The previous section motivated the theory that trion creation or decay drives nuclear polarization effects in our 
system. This theory provides no bias for which direction that polarization is driven, however. In our model, trion 
transitions drive an unbiased random walk in nuclear polarization, with a rate proportional to the probability 
that a trion is created by the pulse sequence. 

Our model therefore is based on a trion-driven diffusion equation. We describe the nuclear hyperfine shift as 
a random variable fl with time-dependent probability density function (pdf) f(fl,t). Time-dependent averages 
over this PDF are given by (g(£l)) — J d^lg(fl)f(n,t), and in particular we write ui(t) = (fi). 

The trion creation probability will be the crucial term in our diffusion equation. We derive this function 
as follows. Let us suppose that, for any pulse sequence, the average value of the electron spin for each cycle 
before optical pumping is (S{Q)) hc{orc , and the average spin afterwards is (S'(ri)) aftcr . Note that coherences are 
not important for the dynamics studied here, so S refers only to S z , and the z is dropped. During the optical 
pumping, the spin is driven to a final equilibrium population S*, as shown in Fig. [2j following a standard Bloch 
equation 

±{S) = -/3(il,\)[{S)-S$, (8) 

where /3(f2, A) is a pumping rate, which depends on the Overhauser shift O and the pumping wavelength A. On 
resonance, this rate has been measured to be on the order of ns _1 in our system!^ After evolving for a time T p , 
we have 

(S(n)W = - e^W) + (S(n)) hcimc e-^ T *. (9) 

This relation may then be combined with a second relation between (<S(f2)) after and (5 , (il)) beforc that depends 
on the pulse sequence. In the single-pulse case, the ir pulse assures (S(^)) a{tcl . = — (^C^)) before' anc ^ therefore 
the one-pulse trion count rate is 

Cx(O) = <S(Q)) after - <S(fi)) before = 2S; tanh[/3(fi, A)T p /2]. (10) 

Note that Ci(f2) depends on A, but we shall not explicitly notate this dependence for brevity. For the two-pulse 
experiment, the pumping laser wavelength is held fixed at A = Ao, but the varying pulse separation of r gives 

<^)>bcfo rc = cosfa, + (6 C + Q)t] (5(n)) after . (11) 

The phase <po depends on offsets in the zero of r. Strictly speaking, at r = 0, the two pulses overlap and 4>o = it. 
In practice, there is a small offset in r that plays little role in the important features of the dynamics. For a more 




Figure 6. (a) The colorscale shows the count rate Ci[u}) as a function of average Overhauser shift u) (vertical axis, 
scale shown in GHz in second panel of subfigure 2), and two-pulse delay r (horizontal axis). The green areas indicate 
where a higher count-rate is expected. Oscillations in the horizontal directions at frequency S a + uj are due to Ramsey 
interference; the Gaussian envelope in the vertical direction is due to the reduction of optical pumping with detuning. The 
superimposed black line indicates stable points where du) /dt = according to Eq. ( 20 1 . (b) Results of the two-pulse model, 



to be compared against Fig. [4] The top panel shows the predicted shape of the Ramsey fringes. These result from an 
Overhauser shift shown on the second panel, which hysteretically follows stable regions of the black curve in (a). The third 
panel shows the pumping efficiency at each r. This particular model uses k/ol = 10 4 ps 2 and /3(w, Ao) = fto exp(— cj 2 /2<t 2 ) 
with a/2% =1.6 GHz and /3 = 3/T p . 



direct comparison to Fig. |4j we set <fio = 0. Hence in the two-pulse case the average number of trions created 
per cycle is now given by 



rmwyrnu maw - z [i - ex P (-/3(n, a )t)][i - cos [(a c + o)r]] 



Now C^fi) depends on r, although again we will not notate this explicitly. This trion creation probability is 
plotted in the green and blue colorscale of Fig. |6ja) as a function of pulse delay r and Overhauser shift O. This 
function oscillates sinusoidally with increasing r due to the spin's Larmor precession, and would lead to fid 
oscillations for fixed SI. The Overhauser shift SI changes the frequency of Larmor precession. The trion creation 
probability drops to for large values of |f2| because the trion transition shifts away from resonance with the 
optical pumping laser, leading to reduced pumping efficiency (/3(f2, Ao) — > 0) and therefore reduced trion creation. 

The three-pulse case, at fixed A = Ao and fixed T, gives a very similar result. The only difference is that 
the relative position of the n pulse drives the rotation angle, rather than the second ir/2 pulse, hence effectively 
doubling t: 

C 3 (fi,T) = C 2 (Q,2r). (13) 



The diffusive terms of our equation are generated by the following argument, supposing a nuclear flip causes 
a change dQ in the Overhauser shift in time dt. The probability of finding an Overhauser shift between Q and 
SI + dfl at time t + dt is, by definition, f(fl,t + dt)dx, and is the sum of three possibilities. In time dt, Jl may 



increase by dfl from fl — dtt, it may decrease by dQ from Q + d£l, or it may stay the same at f2. If Q jumps up 
by dQ with probability [6 + eCj(x)]/2, where S is a trion-independent spin diffusion probability and e quantifies 
the total trion-nuclear spin-flip probability and if f2 jumps down with the same probability, then the probability 
density of seeing f2 at time t + dt is 

f(si ) t+dt) = -[s+ (n + dn)]/(n + dfi, t) + - [J + eCj (n - dn)]/(n - dm, t) + [i - <y - (C 3 (n)]/(n, t) . (14) 

In the limit that (it and dfi become vanishingly small, and defining a = dtt 2 e/2dt and D = dfl 2 S/2dt, this 
expression becomes 



|/M) = ^/(o,*) + aA 



c,-(n)^/(n,t) 



(15) 



The parameter a may be estimated as ^ n r n (^4„/fi,) 2 ss 5 2 /l sec. 

An additional term is required to describe the phenomena by which nuclear spin diffusion reduces large shifts 
and results in a finite- width distribution (of nominal width of order l/TJ .) These effects may be phenomenolog- 
ically captured by a term which drives the mean of /(0,t) to zero at a rate k. The entire evolution for f(x,t) 
might therefore considered to be 



= K~[0/(Q,t)] + ^{[D + 0*7,(0)] \- ( in.) 



Note that if /(fi,i) is a proper PDF with unit integral, then this equation conserves that property. 

Before proceeding, let us note that an equation similar to Eq. ( [16] ) has been used to describe nuclear feedback 
effects in optically controlled QD experiments previously. Reference [7] derives, as Eq. (16) of the Supporting 
Online Material, a diffusion equation similar to the above. However, this experiment, as well as that in Ref. [5J 
are quite different from the one we describe here. That experiment works with an inhomogeneous ensemble 
rather than a single dot, and employs a very different pulse sequence; we have a CW optical pumping field which 
is not present in those other experiments. This pumping field is likely (although not provably) the dominant 
reason the model we present in our manuscript is successful while attempts to apply the models in Refs. [7J[8] have 
failed in our case. A key assumption in Refs. [7j[8]is that the asymmetry in the direction of the average electron 
spin direction during the pulse sequence, or the "light-induced decoherence" of the ground-state electron induced 
by the pulse-modulated electron hyperfine interaction, are key factors in the feedback process. Originally, these 
references led us to believe that this would be true in our experiment as well. However, in our experiment, 
we pump the electron to a particular spin state after each rotation pulse pair, an element not present in the 
experiments described by Refs. 7,8. As a result, the average electron polarization in our experiment is almost 
constant with respect to the parameter r, and therefore does not affect our data, and decoherence of electron spin 
precession is far less important, since only during a very small fraction of time in our one- and two-pulse sequences 
is the electron actually precessing (in contrast to Refs. [7j[8] where the electrons are precessing most of the time). 
We have substantial control over our pulse sequence, and we therefore tested this assumption by varying the 
spacing between pulses, the duration in optical pumping, etc. Doing so allowed us to change the effective electron 
spin orientation during the bulk of the experiment and the relative degree of "light-induced decoherence" and 
we discovered that these factors could not explain our data. In total, our experimental explorations led to the 
following critical conclusion: Only the rate of real trion creation during optical pumping is important for nuclear 
effects in our experiment, independent of the average electron polarization, unlike in Refs. [7j|8] (Of course, in a 
single experiment, the electron polarization right before the optical pumping step and the trion creation rate are 
correlated.) 



We now examine Eq. ( 16 1 in three important regimes. The first regime is where trion-induced effects are 
unimportant. One way to achieve this is to set a <C D; in that case the steady-state solution is a Gaussian 
distribution with mean decaying to zero at rate K and variance which exponentially approaches D/k in time. 
In this respect, we might associate D/k with l/TJ, so that this equation could describe a simple picture of the 
fid; that is, the fluctuation rate would be independent of the pulse-sequence and the fringes of (C(0,t)) would 
undergo Gaussian decay with time-constant T 2 *. Another regime in which simple Gaussian decay is seen is if 



/(f2,0) is Gaussian, and the amount of time t the system is allowed to evolve is small. In this case, solutions 
remain Gaussian and ordinary dephasing processes are observed. This regime is accessed only in experiments 
employing a spin-echo. Critically, in a spin-echo experiment there is substantial time for the nuclei to diffuse 
unhindered by electron polarization, since for the majority of the pulse sequence, the electron spin is precessing 
in its Bloch sphere, which effectively decouples the electron from the nuclear dipole-dipole dynamics. The 
subsequent lack of z-oriented polarization results in a zero average hyperfine gradient on the nuclear spins, 
allowing maximal diffusion and therefore the D 3> a regime. If in addition relatively low averaging times are 
employed (resulting in decay dominated by observable noise), trion effects may be ignored, Gaussian decay may 
be observed, and may be measured, as performed in Ref. 

A second regime is where a ~ D, and a steady-state solution is reached (i — > oo). This is believed to occur 
when diffusion is relatively unhindered, as in a spin-echo experiment, but long averaging times are employed. In 
this regime, this flux-conservative PDE is well-behaved, and predicts that /(f2, t) eventually evolves, from any 
initial condition, into a steady-state solution. That steady-state solution may be written 



f(Ct, t -> oo) = / (0, t -} oo) exp 



o D + aCj(u) 



-du 



(17) 



The prefactor /(0,t — ¥ oo) may be found via normalization. This solution describes a distribution with a mean 
that fluctuates around as r is increased, and whose variance varies periodically with r. Such a solution is 
plotted in Fig. [5^c), upper panel. If this solution is used to find the average expected count rate in a spin-echo 
experiment, (Cs(fT)) as a function of r, the result shows qualitative features observed in spin-echo experiments, 
as seen in Fig. [5jc), lower panel. Oscillations at short times show spiky, non-sinusoidal character, and when the 
oscillations damp at high values of r, they do so to a value close to the maximum count-rate. The physical origin 
of these effects in this model is evidently that, at those places where the fringes look spiky, the nuclei have become 
polarized and their distribution is substantially narrowed. As r is increased, this degree of narrowing decreases 
at the spiky (minimum count rate) locations, but there is little change in the distribution at the maximum count 
rate locations, accounting for the decay to near the the maximum count rate. 

While this solution captures some of the puzzling features of spin-echo fringes, it cannot be a complete model 
for the fid data shown in Fig. [4j or even the hysteresis in Fig. [3j Indeed, a linear PDE such as Eq. ( 16 ) under 
assumptions of full equilibrium cannot show hysteresis. However, only one more simple assumption is needed 
to arrive at a model which does match our data well. This final assumption is that the width of /(fi, r) is 
constrained to be small; much smaller than the inverse Larmor period, despite the influence of trion-induced 
diffusion. This is possible if the nuclear system has neither the time nor the freedom to fully fill the distribution 



predicted by Eq. (16), possibly because nuclear spin diffusion is suppressed by hyperfine magnetic field gradients 



from polarized electron spins. Recall that in our fid experiment, the electron spends the bulk of its time with 
some polarization on the z-axis. This assumption is further motivated by the fact that, in fid experiments, 
the sawtooth-like fringes we observe last for at least an order-of-magnitude longer than T%, suggesting that the 
distribution /(f2, t) is substantially narrowed. 

Under this assumption of a narrow distribution due to suppressed diffusion, we model only the first moment 
of f(Cl,t), oj(t) = (S~J (£)) , and assume the second moment is much smaller than 1/r. Using integration by parts, 



duj 

at 



= —koj + a 



on 



(18) 



Equation (18) is not a closed system, because it still requires knowledge of all moments of /(fi,t) to complete 
the average over Cj(Cl). However, to approximately remove the average, we apply our assumption about the 



width of /(O, t) to write 
leaving, to lowest order, 



Cj[fi] ps C$[w(t)] + [n - w{t)\C'![u{t)\ 

a^ , s 

— = -klo + aCj(u}). 



(19) 
(20) 
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Figure 7. Elements of the model for the one-pulse experiment. The panel on the lower left shows klo (black line) and 
adCi{ui) / dui for several values of A, as A is scanned downwards (indicated by the brown arrows) and then upwards 
(indicated by the blue arrows). The circles indicate locations where kuj = adC\{uj) / dui, i.e. the steady-state solutions of 
Eq. (20 1 for j = 1. As A is scanned downward, the polarization ui "surfs" the same edge of the of the adC\(oS) / 'dui curve, 
and ui steadily decreases, as shown in the lower-right panel. Eventually, at low enough A, the peak of adCi(ui)/dui is 
passed and the only remaining steady-state solution shifts to ui = 0. This solution is maintained as A is scanned upwards, 
until the peak of adCi(u>) /dui reaches w = 0, at which point ui again surfs the adC\(ui) / dui. The resulting countrate, 
C\(ui), is shown in the upper right panel, and is to be qualitatively compared to Fig. [3] This particular model uses 
K/a = 3 x 10 4 ps 2 , 



Invoking our separation of timescales, we presume ui(t) evolves from its initial value (set by the last chosen value 
of A or t) to a quasi-equilibrium final value uif. This final value determines the expected count rate Cj(u>{) at 
this value of A or r. We solve by assuming w(0) = at the first attempted value of A or t, and then we scan the 
parameter as in the experiment, finding the steady-state solution of Eq. (20) at each value. 

The nonlinearity in this ordinary differential equation now allows hysteresis. The hysteresis is most easily 
understood in the one-pulse experiment, whose model is illustrated in Fig. [7] Scanning the parameter A is easy 
to envision graphically, since for this sequence a change in scanning frequency 5ui — S(2ttc/X) amounts to a 
linear shift in C((w), i.e. C[{oj) — > C[(oj + Sui). The steady-state solution is where aC^ui) = kui; several such 
solutions exist at particular values of A. When one such solution lives on the edge of a peak of C[(uj), it "surfs" 
that edge until its crest. This picture explains bistability and hence hysteresis; see Fig. [7| With the exception 
of experimental differences, this model is similar to that used to explain hysteresis in cw experiments!^ 

This basic model allows an analysis of the more complicated two-pulse experiment. The essence of our model 
is that the system quasi-equilibrates to a stable value of nuclear polarization that lives on the edge of the fringes 
shown in Fig.[6ja); with our separation of timescales, the nuclear polarization then "surfs" along the edge of this 
function as r is changed. As r is increased, \u)\ increases causing the observable photon count to decrease due 
to the reduced degree of optical pumping. When \ui\ is so high that pumping is ineffective (/3(to) —> 0) and the 
trion-induced walk stops, spin-diffusion causes the system to drift back to a new stable magnetization at a lower 
value of |cj|, and the process continues. 

Figure [6] shows the modeled Cziwf) and Wf/27r as a function of r. This particular model used ft/a = 10 4 ps 2 , 
which provides the estimate k _1 ~ 100 ms, a relevant timescale for dipolar-induced nuclear spin diffusion. This 
choice of parameters reproduces the qualitative shape of the data quite well, and quantitatively reproduces the 



location where sinusoidal fringes evolve into sawtooth-like fringes. Detailed curve fitting is not reasonable for 
this process, since the quantitative behavior depends on a random initial conditions. Qualitative differences 
are dominated by the random conditions that develop when the stage is moved on its rail, forming the breaks 
between data sets in Fig. |4ja) . 

Despite the small differences, the success of our simple mathematical model in qualitatively reproducing the 
initially mysterious hysteretic and non-sinusoidal phenomena in one-, two-, and three-pulse experiments provides 
high confidence that the basic physics of nuclear feedback is understood in this system. 

6. CONCLUSION 

The effects we explain here, along with similar effects reported elsewhere fDfiEE may be useful for future coherent 
technologies employing QDs. Both the two- and three-pulse sequences may serve as "preparation steps" for a 
qubit to be used in a quantum information processor, as it tunes the qubit to a master oscillator. The efficacy of 
this tuning is evident from the reduced width of /(f2, t) predicted by our model, and appears experimentally as 
spiky fringes in spin-echo experiments and saw-tooth fringes in fid experiments. Tuning to a master oscillator is 
critical for coherent qubit control.^ In the shorter term, the phase- controlled (5-function-like rotation pulses we 
employ are very promising for strong dynamical decoupling (dd)P3 Unfortunately many dd schemes, especially 
those that compensate for pulse errors such as the workhorse Carr-Purcell-Meiboom-Gill (cpmg) sequence, 
require some method to tune the QD's Larmor period to an appropriate division of the pulse-separation time. 
Judicious use of the effects we have modeled here may enable such experiments in self-assembled QDs in the near 
future. 

In summary, we have observed nonlinear nuclear feedback effects in a single charged QD resulting from the 
countering processes of random nuclear walks driven by trion creation, the finite width of optical absorption, 
Overhauser-shifted Larmor precession, and nuclear spin diffusion. This feedback may be employed for tuning the 
electron Larmor period to a particular pulse separation time for more complex pulses sequences. Future work 
might involve extensions of this theory to explain non-Markovian effects of nuclei under more complex pulse 
sequences, as well as exploiting them to extend QD-based quantum memories. 
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